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Using mainly numerical methods, we investigate the 
width of the spin gap of a spin-1/2 two-leg ladder 

described by % = J\Y^ll\ s 'o > ' Sj+i,a + Sj,b ■ 



where Sf 



] = LV\~3,a~j,b ~ u j ,a u j ,b I ~ j ,a u j ,bl> 

denotes the a-component of the spin-1/2 



operator at the j-th site of the a (b) chain. We mainly 
focus on the J r 3> J\ > and |A| <C 1 case. The width 
of the spin gap between the M — and M — 1 sub- 
spaces (M is the total magnetization) as a function of A 



anomalously increases near A = 0; for instance, for 
-0.1 < A < 0.1 when Jy/J, = 0.1. The gap forma- 
tion mechanism is thought to be different for the A < 
and A > cases. Since, in usual cases, the width of the 
gap becomes zero or small at the point where the gap 
formation mechanism changes, the above gap-increasing 
phenomenon in the present case is anomalous. We ex- 
plain the origin of this anomalous phenomenon by use 
of the degenerate perturbation theory. We also draw the 
ground-state phase diagram. 
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1 Introduction Many papers on the spin ladder have 
been published so far 0]. Here we study the ground-state 
(GS) properties of an S — 1/2 two-leg antiferromagnetic 
ladder with Ising-like rung interactions, which has not been 
investigated in detail until now. Our model is sketched in 
FigQ]and expressed by 

N/2 

H = J\ "^^[Sj.a ■ Sj+l >a + Sj.b ■ Sj +i t b] 



N/2 



y][M^ta^tb + Sj a Sj b ) + S z ] a S z j b ], (1) 
0=1 



Figure 1 Sketch of the present model. 

where S" a r b \ denotes the a-component (a = x,y,z) of 
the spin-1/2 operator at the j-th site of the a (b) chain, and 
TV the number of spins which is supposed to be a multi- 
ple of 4. The parameter J\ denotes the magnitude of the 
leg interaction which is isotropic, and J r denotes that of 
the rung interaction with the Ising-like anisotropy which is 
expressed by A. We suppose J\ > and J r > and we 
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restrict ourselves to the J\ <C J r and |A| <C 1 case. Here- 
after we set J r = 1 as the unit energy, In this paper we 
investigate the behavior of the spin gap and determine the 
GS phase diagram on the A — J\ plane mainly by use of the 
numerical methods. 

2 Anomalous behavior of the spin gap Let us 

consider the spin gap defined by Aqi = Eq{AI = 
1,N) - E (M = 0,N), where E (M,N) is the 
lowest energy of the N spin system in the subspace 

M = J2f /2 ( s ta + s lb)- If J i = °. the GS is the direct 
product of the rung state, where the rung state is either the 

singlet-dimer (SD) state [( j) — ( j )] or the triplet-dimer 

(TD) state ^75 [(J) + ( })] depending on whether A > or 
A < 0. The spin gap Aqi is the energy difference between 
the SD state (A > 0) or the TD (A < 0) state and the 
rung ferromagnetic state ( I ), which is equal to (1 + | A|)/2 
having the minimum 1/2 at A = 0. Thus, for sufficiently 
small leg interaction J\ <C 1, it is expected that the spin 
gap Aqi survives and its formation mechanism is different 
for the A > and A < cases. The change of the gap for- 
mation mechanism is often associated with the quantum 
phase transition (QPT). In usual QPT, we see three cases: 
(a) the spin gap Aqi becomes small (but finite) at the first 
order QPT point, (b) becomes zero only at the second QPT 
point, (c) becomes zero in the finite range of the quantum 
parameter for the Berezinskii-Kosterlitz-Thouless QPT. 
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Figure 2 Numerical results of the spin gap Aqi for J\ = 
0.1 with N = 8, 12, • ■ ■ , 28. The N dependence of A 01 is 
reversed between the A ~ 0.0 case and the A ~ 0.1 case, 
although it is very weak. 

Figure [2] shows the numerical results of the spin gap 
Aqi, calculated by the exact diagonalization (ED), for J\ = 
0.1 with N = 8, 12, ■ • • , 28. As can be easily seen, the 
spin gap Aqi always opens and clearly increases in the re- 
gion —0.1 < A < 0.1, which is quite different from the 
above expectations. Thus this behavior of the spin gap Aqi 
is very anomalous. We note that the increase of Aqi in both 
sides are due to the energy difference between the SD (or 
TD) and the ferromagnetic state, which is approximated by 
1/2 + 14 



To explain the anomalous behavior of the spin gap Aqi, 
we use the degenerate perturbation theory considering the 
situation J\ <C 1 and |A| <C 1. First we consider the 
rung state with J\ = 0, which is nothing but the two- 
spin problem. The eigenstates are the ferromagnetic states 
( t) and (|) with the energy 1/4, the SD state with the 
energy — (1 + 2A)/4, and the TD state with the energy 

— (1 — 2A)/4. Here we take the SD and TD states into con- 
sideration, neglecting two ferromagnetic states with higher 
energies. We introduce the pseudo-spin T to express these 
two states. Namely, the eigenstates of T z with the eigen- 
values +1/2 and —1/2, | ff) and | J|), correspond to the 
SD state and the TD state, respectively. Next, by taking the 
effect of J\ into account in the lowest order we obtain the 
effective Hamiltonian 

N/2 N/2 

U eS = 2J 1 J2T*T* +1 -Xj2T*- T (2) 
i=i i=i 
Since both of | ff) and | JJ-) have zero magnetization in the 
original spin picture, TL c s describes only the M = sub- 
space. The 'magnetic field term' (the second term in the rhs 
of eq.©) comes from the energy difference between the 
I ff) and I JJ.) states. The effective Hamiltonian H c s is es- 
sentially the one-dimensional transverse field Ising model 
exactly solved by Pfeuty (2) . The GS of H c s is either the 
disordered state or the Neel ordered state in the x-direction 
(x-Neel state) according as |A| > J\ or |A| < J\. The QPT 
between these two states is of the two-dimensional Ising 
type. In the disordered phase only the expectation value 
(T?) is nonzero; (T/) > for A > and (T/> < for 
A < 0. Thus, the rungs are essentially in the SD state for 
A > J\ and in the TD state for A < — J\. For the .T-Neel 
phase, the eigenstates of 7J are ( j + J]-) ) = ( J ) and 

75H it) + I $)) = (j) with the eigenvalues +1/2 and 

— 1/2, respectively. Thus the a>Neel state in the T picture 
is the z-Neel state of the original S picture. The Neel state 
is realized by the cooperative effect between the SD and 
TD states induced by J\. Three GSs are sketched in Figj3] 
Figure |4] shows the GS phase diagram based on 7i c a- 

The physical interpretation of the above calculation is 
as follows. The energy difference between the SD and TD 
states per rung is |A|. When the perturbation energy J\ is 
smaller than this energy difference, J\ < |A|, either the 
SD state or the TD state becomes important for the rung 
state. On the other hand, when the perturbation energy J\ is 
larger than this energy difference, J\ > | A|, both of the SD 
and TD states are cooperatively concerned with the forma- 
tion of the rung state. The linear combination of these two 
state with equal weight brings about the rung Neel state. 

The Neel state is a doubly degenerate state, while the 
SD and TD states are unique and protected by symmetry. 
Namely, if we exchange the upper and lower legs, the wave 
function of the SD state changes its sign and that of the 
TD state is unchanged. For the Neel state, this operation 
will change a wave function to the other wave function of 
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Figure 3 Sketches of the three GSs, i.e., the SD, TD, and 
Neel states in the original S picture, from top to bottom. 
Ellipses show the SD pairs and rectangles the TD pairs. 
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Figure 5 Behaviors of A$ (left) and Afo (right) for the 

(2) 

J\ = 0.1 case. The quantity A y 00 ' closes in the TV — > oo 
limit at the QTPs within the numerical accuracy. 
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Figure 4 GS of the present model by use of H e s- 



the doubly degenerate GS. The GSs of finite systems are 
usually unique and do not exhibit the spontaneous symme- 
try breaking. For the Neel state, a low-lying excited state 
with M = in the finite system asymptotically degen- 
erate to the GS as N — > oo, and the reconstruction of 
these two states by the linear combination results in the 
doubly degenerate Neel state with spontaneous symmetry 
breaking. Thus, in the Neel region, the lowest excitation 

energy within the M = subspace, Aqq = E\{M = 
0) — Eq(M = 0), will rapidly decrease with the increase 
of the system size N. Here E\ (M = 0) is the energy of 
the first excited state within the M = subspace. The 
left panel of FigfS] shows the numerical results of Affi for 
J\ = 0.1 by the ED, in which we can clearly see the above- 
mentioned behavior for —0.1 < A < 0.1. The true gap 



within the AI = subspace is Aqq in the SD and TD re 

gions, while it is A$ = E 2 (M 
Neel region, where ^(M = 
excited state within the M 



0)-E o (M = 0) in the 
0) is the energy of the second 
subspace. The behavior of 



A 



(2) 

qq is shown in the right panel of Fig. 5. 
Let us define the shift of the GS energy due to J\ by 



AE {M = Q) = E (M = 0,Ji)-E (M = 0,Ji = 0). (3) 
If we use HcH, we can obtain 
AE (M = 0) 
Nil r* /2 

2 I 7T 



2J\\ cos k + X 2 dk 



•(4) 
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Figure 6 Shift of the GS energy AE (M = 0) for the 
J\ = 0.1 case. Dots denote the numerical results with N = 
28 and broken lines the results of the perturbation theory 
eq.©. 



As can be seen in Fig. [6] the decrease of the GS energy 
due to J\ is much larger in the Neel region than in the SD 
and TD regions. The M = 1 lowest state can essentially 
obtained from the GS by changing a rung pair (Neel pair, 
SD pair or TD pair) to the ferromagnetic pair, although the 
location of the ferromagnetic pair is not fixed. Thus, this 
energy cost, which is essentially the spin gap Aqi, is larger 
in the Neel region than in the SD and TD regions. This is 
a physical explanation for the anomalous behavior of the 
spin gap Aqi shown in Fig. [2] We note that we cannot di- 
rectly calculate the spin gap Aqi by use of the effective 
Hamiltonian TL c s, because it describes only the M = 
subspace. 

3 GS phase diagram The boundary between the 
Neel and SD (or TD) phases can be numerically deter- 
mined by use of the phenomenological renormalization 
group (PRG) method |3]. Namely, the critical value A c for 
fixed J\ can be obtained by A c = linijv-^oo X C (N, N + 
4), where X C (N,N + 4) is the numerical solution of 

NA$(N, A) = (N + ±)A$(N + 4, A). Figure shows 



Figure|6]shows the behavior of AE (M — 0) for J\ = 0.1. the P hase diagram of the present model on the \—J\ plane. 

The spin gaps Aqq and A^ within the M = subspace 
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of the present model and the S 
Heisenberg chain model. 



1/2 bond-alternating 



Figure 7 GS phase diagram of the present model deter- 
mined by the ED with the PRG method. Dots denote the 
numerical results and broken lines the results of the pertur- 
bation theory A = J\. 



close on the phase transition lines, whereas the spin gap 
Aqi always opens. Thus, on the phase transition lines, 
at sufficiently low temperatures, the specific heat C(T) 
shows the gapless behavior, while the magnetic suscepti- 
bility Xz {T) for the magnetic field along the z direction 
shows a gapful behavior. In usual cases, both of the spin 
gaps within the M = subspace and that between the 
M = and Al = 1 subspaces close simultaneously on the 
phase transition line, which leads to the gapless behaviors 
of both of C(T) and Xz(T) at low temperatures. 

4 Discussion and concluding remarks Let us 

compare the present model with the S = 1/2 bond- 
alternating Heisenberg chain 



N 

£ 

i=l 



[i + (-iys]s j>a -s j+lta . 



(5) 



where 6 (— 1 < S < 1) denotes the bond-alternation pa- 
rameter and N the number of spins which is supposed to be 
even. This bond-alternation model is a typical model which 
exhibits the change of the gap formation mechanism. When 
S ^ 0, the spin gap Aqi is finite even in the thermodynamic 
limit as A 01 ~ \S\ 2 ^ 3 /y/\\og\6[\ 00. The SD pairs ex- 
ist on the (1 + (5)-bonds for 5 > whereas they exist on 
the (1 - 5)-bonds for S < 0. When S = 0, the GS is the 
Tomonaga-Luttinger (TL) liquid state with the gapless ex- 
citation. Thus the second-order QPT occurs at S = 0. This 
situation is sketched in Fig. [8] As can be seen from Fig. 
[8] the location of the SD pair is different for the <5 > 
and 5 < cases. In other words, the unit cell which is im- 
portant for the gap formation is reconstructed when 5 pass 
through the QPT point 5 = 0. On the other hand, for the 
present model case, the unit cell is always a rung and the 
reconstruction of the unit cell does not occur, as is shown 
in Fig. [3] Only the change of the GS of the unit cell occurs 
when A is swept. This is the very important and essential 
point for the quite different behavior of the spin gap Z\qi 




5>0 

5 = 
TL liquid 

5<0 



Figure 8 GSs of the £ = 1 /2 bond-alternating Heisenberg 
chain Hbn- Ellipses denote the SD pairs. The locations of 
ellipses are different for the 6 > and S < cases. 

As stated above, the key points to the anomalous be- 
havior of the spin gap Aqi are "no reconstruction of the 
unit cell" and "change of the GS of the unit cell". Thus, in 
other models with these properties, similar behavior of the 
spin gap Aqi is expected. In fact, in the 5=1/2 three-leg 
isosceles spin nanotube, a similar behavior of the width of 
the magnetization plateau at 1/3 of the saturation magneti- 
zation is observed (6). 

The anomalous behavior of the spin gap Aqi seems to 
remain even for the the J\ 1 and > J\ 3> — 1 cases. 
We have also succeeded in explaining this anomalous be- 
havior for the former case by use of the bosonization ap- 
proach. The details will be published elsewhere. 

Finally we note that a closely related model was inves- 
tigate by Zheng-Xin Liu et al. \1\ from a somewhat differ- 
ent point of view. 

In conclusion, by numerical methods, we have ob- 
served an anomalous behavior of the spin gap Aqi of the 
spin-1/2 two-leg antiferromagnetic ladder with Ising-like 
rung interactions. Also we have succeeded in explaining 
the physical origin of this anomalous behavior. 
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